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ABSTRACT 


A detailed analysis of temperature variation within 
the sodium-sodium intermediate heat exchanger of the fast 
Breeder test reactor at Kalpakkam has been performed This 
analysis avoids the use of empirically determined heat 
transfer coefficients since it may cause large errors m 
the case of liquid metals The analysis involves a two- 
region Sturm-liouville problem with both negative and positive 
eigenvalues The results of the present analysis predict 
the efficiency of the heat exchanger to be six percent less 
than that predicted by the approximate method using Nusselt 
numbers for constant heat flux* It is shown that the 
approximate method can lead to substantial errors if the 
ratio of heat capacity rates of the two streams is around 
2 (or 0 5) and the heat exchanger is not long 



CHAPTER 1 


INTRO HJCTION 


[11] INTRODUCTION 

Nuclear power competes commercially at present with 
any other conventional source of power like thermal power or 
hydel power There are more than hundred nuclear power plants 
in commercial operation, at present, all over the world 

Nuclear power is produced by controlled fission of 
Nuclear fuel like Uranium- 255 or Plutonium-239 Each fission 
releases energy which is responsible ultimately for generating 
power This chain reaction is carried m a reactor core The 
core is normally composed of fuel, moderator, coolant, control 
devices and structural material Coolant which carries heat 
from hot fuel elements may also act as a moderator m the 
core If moderator is removed from the core and fission 
chain reaction is maintained without slowing down the 
neutrons, then the reactor is called fast reactor Its 
counterpart with moderator is known as thermal reactor 

r 

Transfer of heat to the coolant inside core, is mainly 
by conduction and convection This coolant is pumped to an 







2 


intermediate heat exchanger, (IHX) A secondary fluid picks 
up heat m IHX and goes to steam generator (SG) Steam produced 
m SG runs steam turbine which is coupled to electric genera- 
tor Therefore the major difference in conventional and 
nuclear power plant lies m the mode of generation of primary 
heat 

[1 2] FAST EE ACTORS 

In fast reactors the neutrons are not deliberately 
slowed down by a moderator but only by inelastic collisions 
of neutrons with fuel and other reactor material nuclei, Fig (1 1 
Fast reactors have very high power densities. A typical fast 
reactor has power density about 40 times that of heavy water 
moderated natural uranium thermal reactor 

liquid metals are the only possible coolants for fast 
reactors, because of high power generation per unit volume, 
within the core, liquid metals have excellent heat transfer 
characteristics which results m low hot spot factors, x e , 
less local temperature variations Use of liquid metals 
give high thermal efficiency. 

There is very high neutron flux inside reactor core 
This results m high level of induced radioactivity in the 
reactor materials and liquid metal coolant Coolant is 
circulated m the heat exchanger and therefore there is & 
chance of leakage of radioactive coolant to the other side of 
the heat exchanger 


3 


If on the other side, the working fluid is steam then 
the leakage of radioactive coolant to bhe steam side may 
expose the wo iking personnel (dealing with the turbine 
or condenser) to harmful radiation 

Radio active sodium has a long half life of 15 hrs 
and it emits high energy gamma rays Hence leakage of 
sodium to steam cycle is dangerous This is avoided by 
employing an intermediate heat exchanger (IHX), m which 
liquid sodium flows on both the sides This reduces the 
possibility of radioactive contamination of steam because 
now the steam is generated m steam generator, where the 
secondary sodium (not active) is flowing on the other side 

[12] PAST BREEDER REACTORS (PER) 

Interest m fast reactors stem from the fact that 
they offer a singular possibility of converting fertile 
material to fissile material Past reactors can have ratio 
of conversion of fissile material generated (from fertile 
material) to fissile material consumed, equal or greater 
than one Such a reactor is termed as fast breeder reactor 
( PBR) 

PBR thus offers the possibility of producing power and 
nuclear fuel simultaneously This can be achieved economx- 
cally only when many of the technological hurdles are 


overcome 
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[ 14 ] liquid Mums 

Past reactors do not have moderators and therefore it 
must be ensured that coolant does not act as moderator Hence 
the coolant nucleus must have small cross-section for inter- 
action with neutrons This coolant has to satisfy stringent 
requirements as regards its thermal properties As mentioned 
earlier the rate of heat generation inside fast reactor core 
is very high Therefore coolant must have high thermal 
conductivity, low vapour pressure, high specific heat and 
high volumetric heat capacity These properties will enable 
it to transport heat from core rapidly with less pumping 
power and smaller area of contact Por reactors with high 
thermal fluxes and operating at high temperatures, liquid 
metals have special interest as coolants Liquid metals 
have excellent thermal properties and they fulfil the 
above mentioned requirements quite well 

Amongst various liquid metals sodium is the most 
suitable coolant for a reactor operating at high temperatures, 
particularly for fast reactors where moderating (low mass 
number) materials cannot be used In fact sodium and sodium 
potassium alloys have been successfully used as coolants in 
both fast and thermal reactors. They are also used as 
secondary coolants in intermediate heat exchangers, m 
nuclear reactor power plants 
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[1 5] LIQUID METAL HEAT EXCHANGERS 

Research m liquid metal heat transfer area is going 
on for last three decades Most of it was devoted to develop- 
ing correlations for use m engineering design Tht? Accurate 
and reliable prediction ws-re not considered important till 
recently The mam reason for this was relatively poor 
thermal properties of heat exchanger wall material This 
resulted m the temperature differences across the wall 
thickness an order of magnitude higher than temperature 
differences associated with convection to liquid metals 
Therefore, fairly large inaccuracies (say + 30 percent) of 
liquid metal heat transfer coefficients did not greatly 
affect the overall thermal design But for liquid metal 
heat exchangers used m fast breeder reactors a very accurate 
thermal design is required on account of safety considerations 

Liquid metals have certain peculiar properties In 
liquid metal heat exchangers the use of conventional heat 
transfer coefficient concept is not valid, hence the existing 
methods for nonmetallic fluids are not directly applicable 
to liquid metals 

Liquid metals have small viscosities Therefore 
their application as coolant or heat exchanger fluid, e g, , 
m nuclear reactors, almost always involves large Reynold 
no and consequently turbulent flow, Thermal conductivities 
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of liquid metals are very high, hence molecular conduction 
remains significant throughout the cross-section of flow 
channel or duct geometry Therefore temperature gradients 
are not localized near the wall as it is the case with 
nonmetallic fluids This is a characteristic difference 
between turbulent convection m liquid metals and nonmetals 

Prandtl numbers for liquid metals are very small 
Therefore, Peclet numbers associated with even turbulent 
flow may be as low as 50, with nonmetallic fluids it is 
rarely so Low Peclet number associated with high thermal 
conductivity results m axial heat conduction Millsaps 
and Pohlhausen (1956) and Singh (1957) have shown that upto 
a Peclet number of 100 axial conduction is significant But 
error due to neglect of axial conduction even at a low 
Peclet number may not be large Schneider (1956) determined 
that for Peclet number of 10 and for length to diameter 
ratio of unity, the neglect of axial conduction gave error 
of only 4 percent 

Liquid metal Husselt numbers are very sensitive to 
boundary conditions Sleicher and Tribus (1957) have 
computed the ratio of Russel t numbers for uniform, heat 
flux to that for uniform wall temperature for different 
fluids as a function of Reynolds number This is shown 


NUSSELI, NUMBER RATIO 


FIG 1 2 
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in. Pig (1 2) Prom this figure we conclude that non 
metallic fluids are insensitive to the type of boundary 
conditions specified Metallic fluids on the other hand 
are very sensitive to this and ^variation m Nusselt 
numbers for two different ca e es may be as high as 40 percent 

In liquid metal heat exchangers the heat flux is 
rarely uniform on the other hand boundary conditions which 
specify the wall temperatures are not often accurate 
description of actual situations of practical interest in 
liquid metal heat transfer In certain simple condensers 
or boilers ; wall temperature boundary conditions can serve 
as reasonable approximations, but with liquid metals this 
approximation may be quite inaccurate* Por example, because 
of excellent heat transfer properties of liquid metals , 
conduction through the condenser or boiler walls will be 
much more important than with nonmetallic fluids 

With heat transfer coefficients sensitive to boundary 
conditions and large number of boundaiy conditions that are 
of practical interest, most of the convenience of traditional 
heat transfer coefficient concept is lost Hence one has to 
start from the statement of basic energy equation for the 
given configuration coupled with the careful statement of 
boundary conditions. Assumption of uniform heat flux of 
uniform wall temperature may lead to significant inaccuracies 
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The thermal modelling of I HZ of -FBTE at Kalpakkam 
is considered m this thesis. The temperature profile 
within the heat exchanger is obtained by solving the basic 
differential equations of heat transfer The use of the 
concept of heat transfer coefficient is scrupulously 
avoided This is because it has been shown by £>tem (1966) 
and others that the use of heat transfer coefficient 
concept can lead to great errors m the prediction of tem- 
perature profiles, m liquid metal heat exchangers 



CHAPTER 2 


MODELLING OP INTERMEDIATE HEAT EXCHANGER AT 

KALPAKXAM 


[ 2 . 1 ] introduction 

A fast breeder reactor (PBR) power plant is coming up 
at Kalpakkam, near Madras This plant is expected to 
produce power, commercially, from early mnties Its design 
is based on the Trench PBR called RAPSODIE A detailed 
description of this reactor system is given by Yevick (1966) 
This EBR can be used to breed Plutonium - 239, a fissile 
material, using Uranium - 238 as its parent fertile material 

A prototype of the commercial reactor is likely to come 
into - operation m early eigthies Its designed capacity 
is 45 MW th A detailed description of its intermediate 
heat exchanger is given here 

The prototype is called fast breeder test reactor 
(FBTR) FBTR is a sodium cooled reactor The heat produced 
at the reactor core is removed by circulating sodium through 
the core The primary heated sodium then transfers heat to a 
secondary sodium at IHX«. This secondary sodium then heats 
up the water m steam generator, for running steam turbine to 
produce electric power* 
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There are two IHX m FBTR and they share the load 
equally viz 22 5 MW th Their normal operating condition 
is given m Table 2 1 


Table 2 1 Formal Operating Condition (For single IHX) 


Power exchanged 
Primary fluid 
Primary fluid flow rate 
Primaiy inlet temp 
Primary outlet temp 
Secondary fluid 
Secondary inlet temp 
Secondary outlet temp 
Pressure (inlet to unit) 
Pressure (outlet to unit) 
Secondary flow rate 


22 5 MW th 
Sodium (active) 

125 kg/sec 
5 20°G 
380°C 

Sodium (not active) 

283°C 
514 °C 

5 03 bars abs* Primaiy 

4 81 bars abs fluid 

*74 kg/sec 


[ 2 2] INTERMEDIATE HEAT EXCHANGER 


There are two IHX They consist of a fixed shell, a 
removable portion having the tube bundle and biological 
shield and other provisions for IHX instrumentation and 
control, Fig (2.1) 
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a) Primary Side 

The primary sodium system components have been 
designed m such a way that the main equipments, e g , heat 
exchanger pumps could be removed without disturbing the 
primary circuit Hence the components have a fixed tank 
connected with the primary circuit m which free sodium 
level is maintained with Argon gas cover (See Pig (21)) 

The removable portion is introduced into this from the top 
The fixed shell of the IHX is a stainless steel vessel The 
vessel has a primary sodium inlet nozzle on its vertical 
portion and outlet nozzle on the dished nozzle enclosure 
at the bottom The vessel has double envelop till the 
height upto which sodium is likely to be present Hot 
nitrogen could be circulated through the nozzle arrangement 
around annulus through the shell and its envelop, to preheat 
equipment before charging with sodium 

Electrical resistance type heaters are mounted on 
double envelope of this vessel for heating up this equipment 
during the initial testing of the secondary loop 

b ) Secondary Side 

There are 888 stainless steel tubes of 12/14 mm in the 
tube bundle welded to tube sheets at each end Tubes are 
arranged over 12 pitch circles around a central pipe of 


12 


diameter 224/234 mm ( see Table 2 2) This central pipe is 
also welded to the tube sheets at either end The secondary 
sodium inlet pipe of diameter 200/208 mm passes inside this 
pipe and is welded only to bottom sheet The annulus between 
two pipes is kept under argon atmosphere The differential 
expansion between the central pipe and the tubes and also 
among tubes is accommodated by giving a wave shape (expansion 
bends) to the tubes at the top 

There is a vertical baffle around the tube bundle 
There are inlet and outlet windows for primary sodium flow, 
m this vertical shell The shell is free to move vertically 
during thermal expansion By means of a ring provided around 
this shell, it rests over another ring of fixed shell, thereby 
avoiding by pass of the heat exchanger by primary sodium flow 
There are also a few vertical and circumferential baffles on 
the outer side of the vertical shell arranged around the top 
inlet windows to achieve uniform flow distributions around 
windows and also flow without any gas entrainment problem 


Table 2 2 Arrangement of tubes 


Boy 


of tubes 

Ro-W 

NQ-t .tub££ 

1 


41 

7 

77 

2 


47 

8 

83 

3 


53 

9 

89 

4 


59 

10 

95 

5 


65 

11 

101 

6 


71 

12 

107 

Total No of Tubes 888 , Pitch, 

pitch, circumferential 20 mm 

radial 

19 mm 




13 


There are temporary block devices to prevent incidental move- 
ment of this vertical shell during handling of this equipment 

The entry header for the secondary tube bundle is 
formed by a dished end welded to the top tube sheet Secondary 
sodium enters at the top, flows all the way down through its 
inlet pipe It collects at the bottom dished end before 
passing through tubes It leaves IHX at the top 

The primary sodium enters through the top of the nozzle 
of the fixed shell It gets distributed around m annulus 
It rises above the nnga described earlier and enters the 
shell through windows at the top It flows downwards from 
there around the vertical tubes which carry secondary spdlum 
It finally flows out through the outlet steam nozzle Tlqe top 
dished end collectors and outlet pipe also have electrical 
resistance heater provided for preheating 

c ) Biological Shield Plug 

The plug and the integral part of the removable portion 
of IHX form the closure of this equipment at the top It rests 
over the top flange of this fixed shell and the sealing towards 
the outside atmosphere is assured by means of two '0* rings 
The parts of the heat exchangers which come m contact with, 
sodium are made of Austenitic stainless steel AISI 316 
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[2 5] MODKLlIII(x 

The IHX is a shell and tube type heat exchanger It 
is a single pass counter-flow heat exchanger Secondary 
sodium flows all the way down to the dished end and then 
rises to top exit There is practically no heat transfer 
from or to the secondary fluid during its downward flow 
because of argon shield around it Primary fluid enters 
the annulus around removable shell at a height such that it 
is close to top entrance windows Therefore this heat 
exchanger can be approximately considered to he a single pass 
counter flow, shell and tube type heat exchanger Hot fluid 
(primary active sodium) flowing outside tubes and secondary 
fluid (sodium) flowing inside the tubes 

The geometrical configuration of the heat exchanger 
tubes is such that the problem is one of heat transfer across 
tube bundles for m line flow of fluid This problem has 
been studied by many investigators, employing different 
methods 

a) Equivalent annulus approach was adopted by Dwyer and 
Tu, (I960), Fnedland and Bonilla (1961) and others 

b) Graphical and lumped parameter methods were used by 
Deissler and Taylor (1957), Dwyer (1966) and others. 


HG 2 2 (a) 

Alt dimensions are in mm 



it 
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Characteristic flow region 

I 

* 



FIG 2 2 (b) 
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c) Finite difference method was employed "by Bender and 
Magee (1969) and Oberjohn (1970) 

[2 4J nQUIVALMT AMULUS APPROACH 

In a shell and tube type heat exchanger, the tube 
centres usually make a finite number of equilateral arrays, 
across the heat exchanger cross-section They are individually 
called unit cell or lattice The cell may be a regular 
polygon eg it is rectangular m IHX considered here, 

Fig ( 2 2a) Each unit cell represents typically the thermal 
behaviour of entire heat exchanger 

Once unit cell is recognized, characteristic flow 
region is determined around each tube, Fig (2 2a) It is 
also -of same geometry as its parent cell geometry The 
boundaries of the characteristic flow region are considered 
to be adiabatic If this characteristic flow region is 
circular, the mathematical analysis becomes very much sim- 
plified This is unfortunately not the case with almost all 
shell and tube type heat exchangers In equivalent annulus 
approach this characteristic flow region is replaced by an 
imaginary circular characteristic flow region around the tube, 
Fig (2 2b) Outer boundary of the annulus is considered 
insulated The annular flow area is taken equal to the flow 
area of the characteristic flow region m Fig ( 2 2a) 
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hi jsmg and Eifler (1973) compared results of 
different approaches, mentioned above, along with the 
approximate analytical method they used They have shown 
that when m equivalent annulus approach the ratio of outer 
radius of annulus to outer radius of tube 1 e , p >_ 1 15 , 
the results obtained using equivalent annulus approach do 
not differ significantly from that obtained by any other 
method 

The Kalpakkam IHX has p =1 57, therefore we adopt 
equivalent annulus approach The mathematical simplification 
offered using this method outweigh the insignificant errors 
contained m the final results 




te double pipe countercurrent flow heat exchanger 


CHAPTER 3 


MATHEMATICAL ANALYSIS 

[31] INTRODUCTION 

Employing equivalent annulus approach the shell and 
tube type IHX at Kalpakkam is converted to an equivalent 
IHX where each tube is surrounded by an imaginary annulus 
Outer boundary of the annulus is adiabatic Hence the 
problem to be solved is c ounter -cur rent , double pipe heat 
exchanger where on both sides liquid metal (sodium) is flow- 
ing A line diagram, 3?ig (3 1)» represents typical annulus 
where m regions 1 and 2 secondary and primary sodium are 
flowing Brown (1968) presented mathematical analysis of 
such a system He basically extended the work of Stein (1966)* 
Stem discussed the problem m generality considering symme- 
trical and asymmetrical ducts, co-current and counter-current 
flows, with constant and varying temperature and heat flux 
conditions We will use his procedure to formulate the 
problem Final results are obtained using Brown's modifica- 
tions which simplify the numerical calculations 
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[3 2] ASSUMPTIONS 

The general equation for energy balance for double 
pipe heat exchanger precludes any possibility of solving 
it by analytical means Following assumptions are made to 
simplify the equation 

i) The fluids enter the heat exchanger at uniform 

temperature 

li ) Physical properties are independent of temperature 
m) Frictional heating is negligible 

iv) Axial heat conduction within the heat exchanger 

walls is negligible 

v) Axial heat conduction within the heat exchanger 
fluids is negligible 

•vi ) Velocity distribution of the fluids is independent 
of axial position (l e , fully developed) and is known 
function of radial variable 

vn) The heatisexchanger is m steady state operation 

The first four assumptions are always almost attainable 
m a physical system Fifth assumption for liquid metals 
is valid if Peclet number is greater than 50 The sixth 
assumption is realized m liquid metals because of their 
small viscosities They attain fully developed velocity 
profiles very early compared to thermal fully developed 
profiles For large Reynolds number (Re No >• 10 4 ) the 
velocity profile is nearly flat even close to entrance region 
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[3 3 ] DIMMblONAl PO EMULATION 

The general energy equations for tube and annulus, 
with the stated assumptions can be written as follows 


Tube 




[(1 + 



3T-J r-^z) 


ar 


1 




3Ti( r-^z) 
az 


Annulus 



[(1 



3T ? (r 9 ,z) 

) T> ■ 1 

' 2 ar 0 J 


a. 


•u 2 (r 2 ) ) 


3T 2 (r 2 ,z) 


(3 1) 


(3 2) 


where , 

e H± = Eddy diffusivity of heat in region i 

» Thermal diffusivity of fluid m region l 

l = 1,2 

The above equations are valid both for laminar and turbulent 
flow 

Boundary conditions for these equations are as follows 


1-3 

H 

H 
t— ' 

O 

= T l0 ^ constant ) 

(S 5) 

T 2 (r 2 , L) 

- T 2L (constant) 

(3 4) 



20 


Tube centre line, 


outer onoulus 


3^(0, z) 




= 0 


dTi C l t> 

~~ YU- 


= o 


continuity of heat flux across inner wall 


(3.5) 


r 12 k l 


3l l (r 12’ z) 




r 21 k 2 


3T 2 (r 2 i,z) 


ar. 


(3 6 ) 


conduction of heat across inner wall 


r 12 k l 


^^1 ^ r 12 > Z ^ T a-i 


w 


ar 


1 In ( r 2 l^ r l2^ 


t T i^ r i 2 ' z ^ ~ T 2^ r 2i* 2 ^ 

( 5 - 7 ) 


where , 

k = thermal conductivity of fluid m region i 

l = 1,2 

k w = thermal conductivity of tube metal 

[ 34 ] DIWENSIOKIBSS POEMUIATION 


We introduce following non-dimensional variables, 


A 

r l 

X’i ^ 

- c: — where 0 < r n 

r 12 “ 1 

1 1 

(3 

8) 

A 

r 2 " r 21 

where 0 <, 

v < 1 

(3 

9) 

r 2 

r 22 ~ r 2l 


A 

rjr 

4z 

• where 0 <, 

A A 

z <_ L 

(3 

10) 

z 

Pe l^ 2r 12^ 
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A . A A ^, ( »3)-T 10 A 

T (r ,z) = where 0 < I <1 

21 10 i = 1, 2 (3 11) 

and Pe^ = Peclet number of fluid flowing m region 1 

Substituting these parameters m Eq (3*1) and (3 2), the 
governing energy equations become, 


Tube 


[ f -l ( ^ ) r- 


9r. 


1 


9T l(^l ,z ) 

9*b 


] = gl (? L ) 


Annulus 


— M X 

3Z 


(3 12 ) 


r 2 +a 9 r 2 


7T [f 2 (r 2 )(r 2 +a) 


9T 2 (r 2 , z) 

3r 0 


A 2 9 T 2 ^ r ? ,Z ^ 
" - <'2> “ - 2 - 2 - 


TT 


(3 13) 


and boundary conditions become 
Entrance 

h °) = 0 (3 14) 

T 2 (r 2> L) = 1 (3 15) 

Interior 

afc,(o> z) 

a 


o 


(3 16 ) 
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where, 


i 


K 


3 \( 1 , 2 ) 


ar. 


3T 2 (o, z) 


ar. 


K 


3^(1, Z) 


W 


A 

3r n 


3T 2 (1, z) 


= T„(0, z) - T, (l,z) 


0 


3r. 


(3 16 ) 


(3 18) 


(3 19) 




u 

i 

R 


lav 


H 


¥. 


K 


= 1 + 


-Hi 

a. 


M r i> 


u 


lav 


(3 20 ) 
(3 21 ) 


= Dimensional velocity of fluid m region i 
= Average velocity m region i, 


= 1 , 2 


r 21 //r 22 

R/(l-R) 

c 2^2 

c l^l 


(3 22 ) 
(3 23) 

(3 24) 


S specific heat of fluid flowing in. region i 


= Mass flow rate m region i 


k 


1 il=R I 

^2 11 


(3 25) 
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In 


( ' T 2l / ' T 22^ 


(3 26 ) 


,2 HER 

“ 1+R 


(3 27) 


[3 5 ] SOLUTION OT GOVERNING EQUATIONS 

The Eq (3 12) and (3 13) are solved using separation 
of variables method This leads to the Sturm-Liouville 
problem for two regions involving discontinuity at interface 
Solution of such problem involves positive and negative 
eigenvalues and their corresponding eigenfunctions A 
detailed mathematical analysis is given m Appendix A, for 
the solution of the Eq (3 12) and (3 13) with the stated 
boundary conditions 


[ 36 ] TEMPERATURE DISTRIBUTION 

The general solutions for temperature distribution, are 
as follows 


Tube 


A A A 

R r i » z ) 


O A A 

-B 2 (L-z) 


- °o + I, <R> e n 

n=l 

x 2 A 

+ 0 n B l.a«R) # "“ '] (5 - 28) 


Annulus 


f 2 (r 2 ,z) 


G + 
o 


+ C 


n 


nL lAn Sl ’ n 

h 2,n ( R> e 



, 2 A 
A Z 

n 


n A A 

~P n ( I-®) 

e 


(3 29) 

r 


5 


] 
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where , 



E hA>= 




n-th positive eigenvalue 

Absolute value of the n-th negative 
eigenvalue, 

Normalized eigenfunction m region 1 
associated with n-th positive eigenvalue 
Normalized eigenfunction m region 1 
associated with n-th negative eigenvalue 


G 


n 


Expansion coefficient associated with 
n-th positive eigenvalue 




Expansion coefficient associated with n-th 
negative eigenvalue 


[3,7] El GENFAEUE EQUATIONS 

2 2 

The g eigenvalues X n and ^ are the squares of the 
positive roots of the equations, 


4 (X n ) = 0 


4 (P n ) = o 


(3 30) 
(3 31) 


where , 


=K„ 


da(0, Xj dF( 1, Xj dG(0, X) 


n 


n' 


dr. 


A 

dr-, 


dr. 


— »<!• X n> 


ap(1 - x n ) , 

— 8(0, X n ) 


A 

dr-. 


- IC 


(3 32 ) 
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<rtpj = 


, - K 

n w 


dtfo, P n ) dP*(l, j3j dG(0, pj 


n' 


dr. 


A 

dr. 


+ 


n' 


dr. 


dF*(l» j3 ) 

f*(i, p ) - k r-^- a (o, pj (333) 


dr. 


r n 


1 


The functions P, F*, G*and G- are defined by following 
initial value equations for auxiliary functions 

[38] AUXILIARY FUNCTIONS 


, * A dF ( r, , x„ ) 5 a 

-2- [ f ! ( ^ ) r x ~ — — ] + t ± g 1 (r 1 ) I'(r 1 ,X n ) = 0 


dr. 


dr. 


(5 54 ) 


where , 


dP (0, x n ) 

a- = 0 , F (0, X n ) = 1 


dr. 


f- [f^) r ± ] - p n 2 $ g^) P'*(r 1 ,p n )= 0 

dr l 

(3 55 ) 

= 0 j F*(0, p n ) = 1 

d g (r 2> X n ) 


dr 


dF*( 0 , (3 n ) 


dr. 


[f 2 (* 2 ) ^2 +a ^ 


dr. 


dr. 


^ “ *n w S 2 ( r 2 ^ 


(r 2 + a ) G(r 2 , X R ) 


0 


(3 56 ) 
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dQ(l, x n ) 

dr. 


= 0 , 0(1, X n ) = 1 


dr* * dG* ( r _ , 8) 0 , 

— [f 2 (r 2 ) (r 2 + of) j * 1 ' ° 2 ‘" 2 


dr. 


dr. 


] + ^ « g 2 (r 2 ) 


(r ? + a) 0* (J 2 , (3 n ) = 0 


(3 37) 


d & (1, P„) 

r— ^ = 0 , 0* (1, 6 ) = 1 

dr, n 


[3*9] EIGENFUNCTIONS 

The normalized eigenfunctions are given by the following 
equations, 

dO(o, X J 

h>,n <*L> = C r <U’ VJ/'f'V (5 37) 


B l,n <*L> 


as*(o, p) 

[ r- 3 - F P n )]/Y|<| (3 38) 


dr. 


E, (r 0 ) 
2,n v 2 J 


« u. n ) 

[K S (r 2 , A n )]/f|tf n | (3 39) 


A 

dr. 


E 2,n ^2^ ~ 


dE*(l, P n ) 


dr. 


G^ig, P n )]/Vlir*! (3 40) 


where N n and N* are defined by following equations for 
normalizing factors 
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[3 10] nohmalizing factor^ 


N 


n 


<J0(0 - o ? „ , 

[- -r-®-] B x {[P(r r x n )] 2 ) 


dr. 


p dF( 1 , X ) p „ p 

- HK 2 [ ] 2 B 2 {[S(J 2 , X n )] 2 > (3 41) 


dr. 


and 


K - C 


dG*(0, f3 n ) 


dr. 


V B n {[F*(r lf (3 n )]^} 


P dF*(l, B ) „ A 0 

HK 2 [ ] 2 Bp([G*(r f f3)] 2 >(3 42) 

dr, * * n 


- 2 *• ' 2 f n' 

where n =s 1,2,3» and B 1 £ j- are defined by Bq. ( ^ 35) 
and (A 36) as bulk average temperatures over cross-section. 

[3 11] hXPAN-QION COEFFICIENTS 

The expansion coefficients A^ and C n are defined by 
the following set of infinite linear algebraic equations 


^ B l {B l f n ^ r l^ B l,k^ r l^ } ( 1 " e 


2* 

L 


) \ 


2 

-\ L 


- B 1 lB i,n < r l> B l,k ) °kJ 


= B 1 


(3 43) 
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dl [ - B l (E l.n < J 1> U- 

+ B 1 (E l,n <*!> E l,lc ( ; i) } (I-® 


-^ 2 1 
8 ) \ 
-x 2 £ 
k ><V ] 


= - B 1 ,E l,n <%» 

with n = 1,2,3, 

and for the zeroth order coefficient 


°o = £ l CB 1 {E ^.n <* 1 » 


Ar 


2 A 

- B i { \, k <n» ( 1 - e ' >fc ) -# 

o 


(3 44) 


(3 45) 


where N 0 = l-H 

Numerical procedure to determine expansion coefficients 
is discussed m a later chapter. 


[3.12J GrJdSEBAL HEAT BALANCE 


A heat balance applied to the heat exchanger at any 
point along its length yields following relations for bulk 
temperatures* 

^lav ^ z ^ = H ^ T 2av ^ “ T 2av ( 3 46 ) 

Hence at z a L 

5 xav < B > = H [ J 2av 4) * *2av <°>2 (3 47) 
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Expressing bulk temperatures m series solutions, the above 
equations can be manipulated using orthogonality condition 
Eq (A 37 ) and Eq ( A 34 ) to yield following result 

* laT (a> = (1-H) 0 Q + H J 2av (S) (348 

At the ends of heat exchanger the above equation becomes, 

T 4 (L) = ( 1 -H) 0 + H 

lav v ' o 

and T 2av ( 0 ) = °o 

[3 13 ] HEAT EXCHANGER EFFICIENCY 

The efficiency of counter flow heat exchanger is 
defined as the ratio of actual heat transferred to that 
which would be transferred if the heat exchanger were 


(3 49) 
(3 50) 


infinitely long 

Actual heat transferred 


C 1 W 1 [*10 “ T lav 
° 2 V, 2 


(3 51) 
(3 52) 


Heat transferred for infinite length 


= C i W i ^ T 10 " T 2L^ 


(3 53; 


Ratio of preceding quantities yields 


H < 1 


e = 


T lav< L > 

H 


1 - T 


2av( 0 ) 


(3 54) 


H > 1 


e = *lav(L) = H ^ 1- ^2av(0)^ 


(3 55) 
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or using Eq (3 49) and (3 50) 


H < 1 


P 1 - H r ^ 
e = —— 0 Q + l 


(3 56 ) 


H > 1 


e = (1-H) Co + H 


(3 57) 


physically as length of heat exchanger increases (i e 
approaches infinity) efficiency must approach to unity 

A 

-Hence as L -* ® , e + l Therefore from Eq (3 54) to (3 57) 


A A 


H < 1 as 1 , T lav (L) - H , 


T 2a^°) * °> °o * 0 


(3 58) 


H > 1 as L -*■ ® , Tq av ( H ) -*■ 1 , 
T 2av (°)^ H—l/H , C Q -v 1 


(3 59) 


The series expansions for the end bulk temperatures may be 
written as, 


hav 


= C o + J 1 K B 1 


2 A 
-X L 


T 2av (0 > 


+ B 1 < B l,n < f l” e n ] (3 «» 

3o + 1 [A » B 2 “U ( ^ ,e ' Pn2t 


+ °n B 2 ra 2,n < ? 2>>J 


(3 61) 
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Letting L -*■ « and applying Eq (3 58) to (5 61) “and 
Sq (3 59) to Eq (3 60) , 

00 

H < 1 as 1 ® , l 0 n B 2 {0 2 ^ n (r 2 )> 0 (3 62) 

H > 1 asl+“ j J A n B 1 {E* >n (r^) - 0 (3 63) 

These equations reveal little about individual behaviour of 

the coefficients A^ and C n as L •*> « However, computations 

have revealed that the terms m each of the above summations 
are identical in sign for all values of nj therefore, since 
summations approach zero as L * 


H < 1 

c n - 0 

/s 

as L ->■ » 

(3 64) 

H > 1 

^n •*“ 0 

A 

as L -*■ « 

(3 65) 


[3 14) LONG HEAT EXCHANGER 


For long heat exchangers i e,, L -*■ * Eq (3.64) and 
( 3 . 65 ) reduces series expansions for the temperature distri- 
butions of the heat exchanger fluids 


H < 1 


r l ' 


z) 






-P n 2 (L- Z ) 


(3 66 ) 
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„ * - * -|3 2 (L-z) 

T 2 (r 2 , z) = C Q + B£ >n (r 2 ) e 

(3,67) 
(3 67) 

_ , 2 * 

H >1 I) =C o + ? C n E (i ) e 

n=l 


*2<*2- • °o + J, °n B 2,n e 


(3 68 ) 


\ 2 * 
>A n 2 


n=l 


(3 69) 

For fully developed region only first term of the above 
expansions are required as other terms are negligible m 
magnitude compared to first term 

[3 15] FULLY DiiVLlOPED KEG-ION 

For fully developed region temperature distribution 
is determined accurately enough retaining only first term 
in series expansion 


H < 1 l) - C 0 + A^ ltl (^) e 


-^ 2 (i-z) 


(3 70) 
- & 2 ( L-z ) 

T 2 (r 2 , z) * 0^ + A-,E* n (r ? ) e 


• q 1 2 1 v 2 1 


H > 1 T-^r-L, z) « C Q + (i^) e 


2 * 

“ X 1 2 


(3 71) 


(3*72) 


2- 

* , " X 1 2 


®2^ 1> 2 > z ) s G 0 + C l Jj 2,l ^ r 2 ^ e 
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Liquid metals have large thermal entry lengths compared 
to nonmetallic fluids Fully developed region for liquid 
metals may not he realized for small or moderately long 
heat exchangers Even for long heat exchangers, for accu- 
rate description of temperature at entry and exit, more 
number of terms are required 


CHAPTER 4 


PjjUG PLOW idealization 

[ 41 ] INTRODUCTION 

In principle assuming that sufficiently accurate 
values of f (r ) and gj ( r ) can be found from recent experi— 
mental results for liquid metals, the Lq. (3 34) to (3 37) 
can be integrated with stated boundary conditions Thus 
auxiliary functions P, P*, G and G* may be obtained This 
may not be straight forward as the differential equations 
may not submit to analytical solution. Expansion coefficients 
will have to be integrated nume rib ally Obviously this would 
not desirable if it could be avoided* Much of engineering 
research can be associated with attempts to discover simpli- 
fications and convenient computational procedures 

Por turbulent flow of fluids, considerable simplifi- 
cations are offered if plug flow model is assumed Then it 
does not become necessary to resort to numerical techniques 
for determining auxiliary functions, m our problem 

[4,2] PLUG PLOW MODEL 

The plug flow model involves the assumption of uniform 
velocity profile i e , ^ 1 and if m addition we neglect 

effect pf eddy diffusivity l e , f^r^) = 1 > ^ 13 possible tp 



RELATIVE THERMAL RESISTANCE 


FIG 4 1 



RELAH/E DISTANCE FROM WALL 


Relative 1 mal resistance 

' >lds 
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obtain analytical solution of the Eq (5 34) to Eq (3 37) 

With these approximations the equations for energy balance 
take simplified form for region 1 and 2 In fact with plug 
flow and f 1 (r 1 ) = 1, Eq (3 34 ) to Sq (3 37) take form of 
familiar Bessel equations and a relatively simple expression 
is obtained for F, F , G- and G* 

Heat transfer analysis made with plug flow model for 
nonmetallic fluids give inaccurate results Stem (1966) has 
shown that with liquid metals this modelling yields meaningful 
and fairly accurate results, for low Peclet number (Pe 4 . 100) 
Peolet number of 100 normally corresponds to a Reynolds number 
of 30,000 This range of Reynolds number for heat exchangers is 
definitely of practical interest Plug flow idealization is 
appropriate for liquid metals because steep temperature 
gradients are not localized near wall so that the error 
resulting from using uniform velocity profile 1 e , g 1 (r l ) = 1, 
in this region is not significant For low Peclet number 
(Pe 4 100) the assumption of f 1 (? 1 ) = 1 1S appropriate because 
eddy conduction plays a minor role compared to molecular 
conduction m the transfer of heat This is demonstrated 
by Fig (41) In our case Pe > 100 and hence an extension 
of this model will be needed This is discussed m Sec- 
tion [4 11] 
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[4 3] PLUG plow solutions 

With = f 1 (r l ) = 1, Eq_ (3 34) to (3 37) become 


_a_ (- "<*!> *n> 
1 aJ^ 


) + h » ( ; i> - 0 < 4 15 


ar(°, \) 

— = 0 , p (o, x n ) = 1 


A 

dr-, 


4- (J, 

dr i 


as*(r 1( B ) 


ar. 


— ) - P n 2 h 0 n ) = 0 < 4 2 ) 


M*(o, P n ) 


dr-, 


= 0 , P* (0, p n ) = 1 


“7“ ( (r 2 + a) 
ar 2 


dG(r ?> x n ) 

■ i fa ’ » » ^ 

dr„ 


- X n 2 “ 2 (V° > e(J 2> »n> = 0 


dG (1, X ) 

-a- = o i g (i, ^) = i 


A 

dr. 


, dG*( r 0 , p„) 2 2 

d / t * c 'ns, o t c- 


dr. 


( (r« + ct) 


dr. 


) + P n *> ( r2 + o ) G (r 2 , P n ) = 0 

(4 4) 


d0*(l, P n ) 


dr. 


0 f 


G*(l. Pq) 
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Solutions of Eq (4 1) to (4 4) m terms of F, F* , G 
and G* are, 

S' (r x , x n ) = J 0 ( X n r x ) (4 5) 

& = J 0 d) (4 6 ) 

G ( r 2* = ^ X n ( 1+a ) (K^[ui X n (l+ a ) ]I 0 [a> X n (r 2 + 0 ^ 

+ IiCw^l+a)] K 0 [u>x n (r 2 + cr)]) (4 7) 

^(r 2 . %) = -f^ n (1+cr) (J 1 [«^(l+a)3T 0 [«|3 a (r 2 +cr ) 3 

+ Y^o^d+cr)] J 0 [“p n (r 2 +a )]) (4$) 

where , 

J l* V Y 0 Bessel functions of first and 

second kind 

If, I Q , Kj, K Q Modified Bessel functions of first 

and second kind 

Substitution of above expressions into Eq (3 30) and (3.31) 
yields eigenvalues Eigen functions can be obtained from 
Eq. (3 37) to (3 40) Finally expansion coefficients from 
( 3 43) and (3 44) 
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[44] NARROW ANNULUS APPROXIMATION 

A useful simplification is obtained letting R -*■ 1 This 
physically means & narrow annulus. As R + 1 , c + » We 
find that with E = 1, F and F* remain unaffected Bq (4*7) 
and (4*8) can be simplified letting a •* « and retaining only 
first term m large arguments approximations for Vessel func- 
tions A simpler approach is letting a ■* » in Bq (4 3) 
and (4 4) which, leads to following equations 


d 2 G(r 2 , A n ) 


dr. 


a. 2 x n 2 s(r 2> X n ) 


0 


(4*9) 


a2 &< G 2 , p k > t 

dr 2 



with initial conditions remaining unchanged 


0 


(4 10) 


It was found thet by both techniques the expressions 
for functions G and G* obtained were exactly same Numerical 
computations have revealed that for R as low as 0 5, error 
m computing efficieppy is less than 3*5 percent For 
present case R « 0 6^6, therefore we adopt narrow annulus 
approximation approach The simplifications offered by 
employing narrow annulus approach can definitely be used to 
solve the problem, because the errors contained m the final 
results are insignificant 
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[4 5] AUXILIARY RUNG HO US 

For narrow annulus approximation the auxiliary 
functions obtained are as follows 


1 'V 

i n ) 

* J o < x n 

(4 

11) 

F *(n- 


- h <Pn> 

(4 

12) 

® (? 2 > 

*n> 

= Gosh [wX n (1 ~ r 2 )] 

(4 

13) 

®*(r 2> 


= Cos OP n (l-r 2 )] 

(4 

14) 


For narrow annulus approximation the auxiliary functions are 
considerably simplified. Also R is eliminated except that 
it occurs m nondimensional parameter K 

[4 6] eigenvalue equations 

Substituting F, F* , G and G* in eigenvalue E a (3 30) 
and (3 31 ) we obtaih, 

4 (X n ) = « \ Slnh ( “ X n> '“ l n J o (l n ) 

Bxnh («X n ) + K Vl (> n> oosh (4 14> 

and 

4 *( p n ) = .K* P n Z I x (P n ) Bin + “ p n V'R 

sin U^) - Ii (%) 005 ( “^n 5 (4>15) 
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[4 7] NORMALIZING PAG TORS 

From equations (3 41) and (3 42) normalizing facto rs 
are given as, 


N. 


n 


•V sinh 2 ( w x n ) [J 0 2 (^ n ) + J;L 2 U n )] 


1 ? 5 Sinh (2 <I>X ) 

-i IiK 2 X n 2 J x (x n ) [1 +-2771—^ " ] 


n 


(4 16) 


N n 


u) 2 p n 2 sin 2 (w^) [I Q 2 ( (3 n ) - I 1 2 (p n )] 


sin (2 »8 ) 

? p n % < t 1 + ~ tv r*~ ] 


1 TTT/-2 n 2 


n 


(4 17) 


[4 8] NORMALIZED EIGENFUNCTIONS 


E i, n ( n> = 

- w x n Sinh (wx n ) J 0 (x n r^/VlNj 

(4.18) 

= 

" P n ain l o <fJ n ? i>/VlK‘l 

(4*19) 

E 2,* n = 

" Kx n J l^ x n^ C08h C“ l N J 

(4* 20) 

B 2,n = 

K ^n cos U-^)] /Y 

(4 21) 


[4.9] EXPANSION COEFFICIENTS 

In Eq (3 43) and {3 44) various factors for plug flow 


are given as follows 
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B 1 {B l,n (r l> B ^, k ( f i>> = 2/ (r x ) \ to. 


.2 0 2 __ „2 / .. „ \ 2/ 


P n Sln <“^n> IX W - *1 X>] 


2__LLS_ 




2 A 

+ 8 L 
, K n 


In* 


n 1 


+ fti^D 


1-e 

for n ss k (4 22) 


2 S xn(o,p n )][ 

- PiXlX 5 KX 


for n ^ k ( 4 23 ) 


B i a i,n (? 1> Sl.k^l)! = 2 / L i,n B l,k < ? 1> a *i 


~ - 2 Oj3 n sin (wp n )] [ “* k sinh ( uX k )][P n I i^ n ) J 0 ^ X k^ 

+ \ I 0 (?„> WtfltVQ l H kl > <Pn 2 + X k 2)) (4 24 > 

B 1 tE l,k < r l> E l,n = 2 i E l,n (r 4 E l,k (r l' r l ir l 


- 2 [ojp k sin (uj^)] [wX n sinh (»X n )] [ 1^ ( J 0 C * a ) 

+ x n i 0 c^) v^i/twy i»ki >x 2 + (4 - 25) 
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B 


X {B l,n E i,k <^i> *> 2 / B 1>n (f x ) E 1(k dr 


C“ x n Slnh <“ x n ):]2 + J i 2 ( » n > ] / |Hj 


+ e 


2 a 

-XL , .. . 

n /(l- e n ) 


-A„ 2 L 


for n = k 


(4 26) 


* 2 ^ x n sinh <» x n )] ^ x k sinh (wX k^ t x n j l (x n ) W 


- x k J B x k> J 0 ( x n )3/C(VlHj |H k ! )(x k 2 - * n 2 )] 


for n ^ k 


(4 27) 


Bi {Biia a x h 


- 2 / B l,n <*L> ? 1 


= 2 


M g n sin( M p n ) I x (p n ) 


(4 28) 


(ri<i ) fi n 


B X {:B l,n < ? 1» = 2 / B l,n < ? 1> *1 


sxnh Ux ) J, (x w ) 


= - 2 


n 


- / w *| V A v- / 

£ 1 — £~ ( 4 29) 


<*!»„! > v 
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mal resistance distribution 


45 


[4 10] BULK TEMPERATURE 

Expressions for bulk temperatures are given as, 


T lav U) 


wem (u>8 )l,(S) -B 2 (L-z) 


= C n + l C2 A. 


n 1 p n 


n 


n=l 


n 


V|H*| 


-20 


wsinh (ui\ ) <L ( A ) -X z 


n 


Vlisr. 


V-, / \ "y, / 

- e n ] (4 50) 


n 


T 2av - 


C o + l c^n 

n=l 


2Ti(p n ) sin(wp n ) -p n 2 (l-z) 


“ f|N* 


K sinh (uX ) J, ( A ) -A 




n 


U) 


r |n. 


n ] (4 31) 


n' 


It can be verified that the orthogonality condition is 
satisfied for plug flow, narrow annulus (R = 1) case General 
heat balance is also satisfied 

[4 11] plug PLOW EXTENSION FOR LARGE PECLET NUMBER RaJJGE 

Liquid motal heat exchanger applications to nuclear 
tecnology involves very high flow rates Conseouently 
Reynolds numbers are large (Re No £.10^) Peclet number 
corresponding to such range is much greater than 100, 

For this range of peclet number, m plug flow idealiza- 
tion, assuming f x (r ) =1 is inaccurate For Reynolds number 
of 10 6 there is a large variation in eddy diffusivity across 
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duct cross-section as is seen in Pig (4 2) Stem (1966) 
suggested an approximation to extend the validity of plug 
flow for peclet number range upto 1000 This approximation 
involves the use of a properly averaged space independent 
total conductivity which results when a simple variable 
change is made 


The following derivations follows the method suggested 
by btem (1966) 

[4 12] DERIVATION 


A variable change m radial variable r is introduced 
The new variable x (r ) is related to r 1 by a first order 
differential equation 


Tube 




, <r+ . 


(4 52) 


Annulus 



x 2 + a 
r 2 + a 


+ 



f 2 (r 2 ) 


(4 33) 


with the boundary condition x x (0) = 0 and with k^ a constant, 
so chosen that x x (l) = 1 
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Introducing variable change m the equations (3 12) 

and (3 13)} m terms of x ±f they can now be written as 
follows 


lube 


1 3 





( i> 


3l^( X^ , Z ) 

. *• 

3z 


(4 34) 


Annulus 


(x 2 + o) ax 2 


3 ... , 3T J ( x 2 . B) 1 

[(x 2 + o) — £ — ] 


3X 2 


.2*1 


- O > 

- «' g 2 (x z ) -M 


3T 0 (x o ,z) 
. < 

■M 

ax, 


whe re , 


z a k| z ; and L as 1 


A A 

r. p dr n 

s i( x i) K g i^i^ ^ 7T 

x i ax i 


r 2 + a N 2 dr 2 


g P (x 2 ) = g (r 2 ) (~ ) 

* 2 ^ x„ + a 


dx. 


(4.35) 


(4 36) 
(4 37) 


(4 38) 


The boundary conditions become, 
Entrance . 


T 1 ( V 0) 


0 


(4.39) 
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<r* iM •— 

T 2 (x 2 , L) = 1 

Interior 

3 ^1 ( ® * z ) 

= 0 

3x l 


(4 40) 


(4 41) 


K 



8T-L (1, z) 



3^2 (°> z ) 

3X 2 


(l , z) 

\ k l = = T (0,z) - T (l,z) 

8X 1 x 


(l,z) 

— 4 = o 

ai 2 


(4 42) 


(4 43) 


(4 44) 


For turbulent flow and all ranges of Pe No the g(x^) and 
g(x 2 ) as defined by Eq (4 37) and (4 38) have average value 
over cross-section equal to unity For small Pe No g(x^) 
and g(x 2 ) are approximately equal to unity for a large portion 
of cross-section of the heat exchanger duct 


Liquid metals are less sensitive to the exact represen- 
tation of velocity profile compared to eddy diffusivity 
representation, Stem (1966) Therefore as an approximation 
it will be assumed that g(x- L ) and g(x 2 ) are identically equal 
to unity for all values of Peclet No. This reduces Eq (4 34) 
and (4 35) to following 
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Tube 


t 





(x lt z) 
ai l 


Annulus 


z) 

3z 


(4 45) 


V[(J 2 + 

X 2 + or 3 X 2 


3T.(x , z) 

a) — ^ 

3x 2 



8T 2 (x 2 , z) 


3X 


2 

(4 46) 


with boundary conditions remaining unchanged 

Eq (4 45) and (4.46) are lusntieal to Eq (3 12) and 
(3 13) when plug flow model 1 e , f 1 (r l ) = g-Jr^) « 1 is applied 
to later set of equations, and making following changes m 
non-dimensional parameters 


kt 

K - — K 

k 2 

(4 47) 

^w '* ^w 

(4 48) 

^ 1 A 

z .+ k^ z 

(4 49) 

In effect, the k* approximation has changed molecular conduc- 
tivity k i to k*, in plug flow equations Hence the results 

obtained for plug flow without k* are valid for large Pe No 


approximation with above modification m parameters only 
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[413] Determination of k* 


eqns 


Expressions for k* may be determined from integrating 
(4 32) and (4 33) between 0 to 1 




1 

J 

0 

*1 (r x ) 

1 

e 

1 ? »1 

i 

^1 ( r i ) 

i 

dx 2 (r 2 ) 

J 

(x 2 (r 2 )+ a) 

i 

( 

<u 

i 

i 

0 

i ..... . ^ ... 

f 2 (r 2 ) (r 2 +o 


(4 50 ) 


(4 51) 


There are two complications m determining k* First it is 

necessary to know f^rj which essentially is knowledge of 

e as a function of radial location This can be obtained 
Hi 

from experimental results Second complication is that right 
hand sides of equations (4 50) and (4 51) take indeterminate 
form at r^ = 0, i e , x x = 0 

An alternative method is suggested by Stem (1966), 
according to which liquid metal he a t transfer coefficients 
can he computed from plug flow values by simple relation 


N u (fc) = k + [N u (z)]pp 


(4-52) 
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where, 

PP - Plug Plow 
Nu - Nusselt Number 

Since detailed temperature distributions are not required 
there is no need to determine x^r.^) f rom equations (4 32) 
and (4.33) Appropriate values of k^ + can be determined 
from equation (4 52) when applied to fully developed 
nusselt numbers either from experiment or by computation 

The equations used for evaluating k.J*' are , 

Tube 

k x + * (7 + 0 025 Pe° 8 )/8 0 (4 53) 

Annulus 

k + =* (5 25 + 0 D188 Pe 8 8 (^)° ^)/6 0 (4 54) 

* r 21 

Eq (4 53) is obtained from correlations given by Lyon (1955) 
and Eq (4.54) from correlations given by Bailey (1950) 


CHAPTER 5 


RECULTS AND DISCUSSION 


[5 l] EIGENVALUES AND EIGENFUNCTIONS 

o p 

The square roots of the eigenvalues, A n and P n , 
were determined from the roots of the Eq (4 14) and (4 15) 
These nonlinear algebraic equations involve nondimensional 
parameters R, and The values of these and other 

parameters for the Kalpakkam heat exchanger are given helow 
To calculate these values thermal conductivities of sodium 
and austenitic stainless-steel were obtained from Rasnjewi 
(1976), 


Pe 1 

car 

160 

Pe 2 

as 

190 

k l + 

= 

1 055 

k + 
2 

a 

1.11 

K 

= 

( 5722) (k 1 + /h 2 + ) 

H 

ss 

1.65 


*s 

( 5698) (k 1 + ) 

% 

ss 

(2 6 )(k^) 

The values 

of 

kj* and k* were obtained 


and (4,54), Numerical values given within 


from Eq (4 53)* 
the brackets are 


plug flow values 
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The Eq (4 14) and (4 15) having &, and 00 as 
multiplying constants were solved by Regula-False method 
The accuracy m determining the roots of these equations 
is of the order of 10 Once the eigenvalues are known 
it is easy bo calculate eigenfunctions from Eq (4 18) to 
(4 21) 

[5 2 ] EXPANSION COEPPICIENTb 

The expansion coefficients were calculated from the 
Eq (4 22) to (4 27) and the nonhomogeneous part from 
Eq. (4 28) and (4 29) To obtain and an infinite set 
of equations must be solved It has been shown by Brown (1968) 
that it is enough to consider only ten equations to obtain 
C 0 accurate up to three decimal places 

Bor long heat exchangers 1 e , z ^ 1, only first term 
was required to calculate temperature distribution The subse- 
quent terms had negligible contribution except close to the 
entrance or exit of the heat exchanger 
[5 3] COMPARISON OP PRESENT METHOD WITH NTU-MEIHOD 

The performance of a conventional heat exchanger is 
usually predicted by the effectiveness - NTU method In this 
method, the calculation of NTU (Number of heat transfer 
units) involves the estimation of heat transfer coefficient 
for primary and secondary side, using ns 

~~r,rwff> 1 

ACC. No £ VS 4#.,<s<ss. 
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As already pointed out in Chapter 1, the heat transfer 
coefficient is a strong function of boundary conditions 
There are two extreme boundary conditions normally considered 
in heat transfer problems, viz , constant heat flux and 
constant temperature 

In the case of flow inside tubes experimental corre- 
lations are available for both conditions but for flow m 
annulus, correlation for constant heat flux is only avai- 
lable. Hence we compare our results with the results obtained 
by NTU-Method, for constant heat flux boundary condition 
Table 5 1 compares the efficiency and exit temperature values 
for primary and secondary sodium as obtained by present method 
and uniform heat flux method For uniform heat flux (UHF) 
calculations see Appendix B, 

Table 5 1 


» 



Present 

Method 

Manuf acturer* s 
Data 

UHF 

Method 

Efficiency 

95 % 

97 5 % 

99 3 % 

Exit Temperature for 
Primary Sodium 

383°C 

380°C 

377°G 

Exit Temperature for 

Secondary Sodium 

508°C 

514°C 

5t8°0 
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The UHF approach over predicts efficiency, but the 
error is not serious because the actual efficiency quoted 
by manufacturer is close to hundred percent It is also not: 
very clear that manufacturer's data represents experimental 
values or calculated values Table 5 1 also shows that UHF 
over predicts exit temperature for secondary sodium and 
under predicts exit temperature for primary sodium The present 
method shows the opposite trend Temperature variation for 
UHF and present method along length is shown m Fig (5 1) and 
Pig (5 2), The maximum temperature difference m two 
methods are at the outlet since the inlet temperatures were 
assumed to be same m both cases 

Manufacturer's quoted temperature values are only 
available at present Actual temperature values will be 
known only after experiments conducted at Kalpakkam These 
experiments might throw light upon the actual obtainable exit 
temperature values At present the Intermediate heat exchanger 
is under fabrication by BHEL 
[5 4J EFFECT OF VARIATION OF FLOW RATE 

Under noimal operating conditions for IHX, H » 1 65* 

A change in operating reactor power level will have a corres- 
ponding change m the quantity of heat transferred at IHX 
from primary to secondary sodium. If reactor power 

level goes up more heat will be exchanged at IHX Consequently , 
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a higher secondary mass flow rate will "be required to 
transport this heat If reactor power level goes down 
less heat is exchanged and a lower secondary mass flow 
rate will be required In either case H will vary as it 
represents ratio of primary sodium to secondary sodium 
heat capacity flow rates 

It is reasonable to assume that primary flow rate and 
its exit temperature from the core is constant When power 
level goes up, the higher power is transferred to primary 
sodium inside core This can be achieved keeping same 
mass flow rate and exit temperature, from reactor core, with 
a lower inlet bemperature to the core When power level 
reduces, other conditions remaining constant, a higher inlet 
temperature to the core is anticipated 

Secondary sodium mass flow rate will depend upon the 
quantity of heat to be transferred But it is reasonable 
to consider exit temperature of secondary sodium from steam 
generator to be constant This is possible because a 
changed secondary mass flow rate will result in a change m 
the quantity of steam generated for same exit temperature 
from steam generator Only IHX exit temperature will 
increase or decrease depending upon increase or decrease m 
reactor power level This may be of practical interest also 
if m secondary loop due to blockage or additional resistances, 


FIG 53 



VARIATION OF PRIMARY SODIUM 
OUTLET TEMPERATURE vs RATIO C 
MEAT CAPACITY FLOW RATE , H', 


TEMPERATURE (°C) 
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the secondary sodium mass flow rate changes Any change m 
secondary sodium flow rate due to either reason, will result 
m a corresponding change m H With the reasonable assump- 
tions of keeping inlet temperatures to the IHX constant, 
exit temperatures for both primary and secondary sodium 
can be predicted, under steady state conditions Fig (5 5) 
and Pig (5 4) gives such a plot, both for present method 
and UHF method 

As anticipated for higher values of H both predict 
higher exit temperatures for secondary sodium For H decreas- 
ing both methods predict decrease m secondary sodium exit 
temperature This trend is observed for primary sodium also 

It is interesting to see that for H > 1 UHF over predic 
the exit temperature values for secondary sodium and under 
predicts for primary sodium, which essentially means over 
prediction of efficiency compared to the present method For 
H < 1, it under predicts efficiency A good agreement 
is obtained for H s 1 This means that for H = 1, UHF method 
and present method yield same results 

For H > 1 and increasing H will decrease C m:Ln , which 
for H > 1 is given by secondary heat capacity flow rate 
Since NTU is defined as W c mn » decreasing C min will 
increase HTU (See Appendix B) and from Eq (B 6) it is obvioi 
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that efficiency will increase This increase m efficiency 
is much higher for UHF compared to present method for equal 
increase in H. Hence UHF method over predicts temperatures 

For H < 1 and decreasing H will fix C „ because now 

min 

^mm 00 ^responds to primary heat capacity flow rate which 
is considered constant From Eq (B 5)» it is clear that 
denominator does not decrease as fast as numerator Conse- 
quently efficiency increase for H < 1 and decreasing H is 
at a slower rate than for H > 1 and increasing H This 
resulted m under-prediction of efficiency for H < 1 by UHF 
method compared to present method 

For Ii * 1, for a heat exchanger energy balance equa- 
tion at any point along the length predicts constant bulk 
temperature difference between primary and secondary 
sodium This means that there is a constant heat flux all 
along the length Therefore UHF should predict accurate 
results under such condition Present method also yields 
same results and therefore it is advisable to use NTU method 
when Hal, based on UHF conditions 

[5 5J COMPARISON OF k + WITH PLUG- FLOW MO EEL 

The Fig. (5 1) to Fig (5.4) were presented using k^* 
approximation. Neglecting k x + in the present case results 
in the under prediction of primary sodium exit temperature by 


i- 



FIG 5 5 



MANUFACTURER S’ DATA 


TEMPERATURE VARIATION ALONG 
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3 0 a nd over prediction of secondary exit temperature sodium 
ty 4 C The Fig (5 5) shows predicted temperature profiles 
along heat exchanger length for the two cases, viz , with k + 
and without approximation The peclet numbers for the 
present heat exchanger are not very for removed from 100 
Therefore, the effect of eddy conduction incorporated using 
^1 n0 ^ vszy pronounced The two temperature profiles 

are therefore very close to each other except at the ends 
where maximum difference is observed 

[5 6] VMI ATI Oh OF HEAT FlUX AL04G- HEAT LXCHMG&K LEHGrlH, 

The heat flux for the present heat exchanger varies 
strongly along length From Flg t (5*6) we see that the 
heat flux at one end is ten times at the other end Hence it 
is incorrect to assume uniform heat flux boundary condition 
for present heat exchanger But results obtained for the press 
case even assuming UHF do not lead to serious errors This is 
because for H > 1 the over-predicted value of efficiency is 
near unity and quoted value is also very close to unity 

For short heat exchangers with smaller efficiency and 
H ^1, UHF assumption may lead to gross over prediction or 
under prediction of the heat exchanger performance But when 
H * 1, the present method may yield results which may be very 
close to UHF predictions 
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[5 7] NAHA! TMPAlUIUfifi VARIATION 

UHF method can predict only bulk temperature variatio 
of the heat exchanger along the length But present method 
was evolved solving fundament al equations for energy balance 
The solution thus gives temperature at all points within the 
considered boundaries Hence present method provides more 
information than UHF method 

A detailed knowledge of temperature at all points, 
is of interest to the designer The radial temperature 
profile gives the magnitude and location of hot spots The 
Fig (5 7) shows that at the centre of heat exchanger, 
maximum temperature occurs at the outer boundary of the 
annulus The present method can be used to determine such 
radial temperature profiles, at any distances from either 
ends* 

Only first term is required to calculate the radial 
temperature profile m fully developed region, i e , m the 
middle of heat exchanger*. But close to the entrance and 
exit more number of terms are required Number of terms 
depend upon the accuracy with which temperature is to be 
predicted Foi short heat exchangers fully developed region 
may not be realized at all This requires large number 
of terms to calculate radial variation along the length 
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Prom Pig (5 7 ), it is found that the temperature drop 
m tube metal is maximum, i e , 17°C Temperature drops m 
primary sodium is 11°C and m secondary sodium it is 16 5°C 
Temperature drop is proportional to the resistance offered 
by the medium Prom Pig (5 7), we conclude that for the 
present case resistances in all three regions are of the 
same order of magnitude 

[5.8] SCTGGLSTlijD 1 UTUHS WORK 

a) Recent eddy diffusivity and velocity distribution 
correlations for liquid metals can be made use of m deter- 

•j* 

mining k 1 accurately It was not required for the present 
problem because the Peclect numbers were not much greater 
than 100 

b) The present method involves use of plug flow model 

where we assumed f 1 (r ) = = 1 This assumption 

can be removed if experimental correlations available for 
flow along rod bundles and flow inside tube are employed to 
determine f^r ) and g 1 (r 1 ) The basic energy equation 

can then be solved numerically as a two point boundary value 
problem Such a general solution will be of immense value 
It can be also used to ascertain validity of experimental 
correlations and validity of any approximate theoretical 
approach like equivalent annulus approach and RTG method 


60 


BIBLIO GRAPH! 


♦ 

t 


i 

$ 



1 

2 

3 

4 


5 


6 


7 

8 


9 


10 


11 


Bailey, R V , 1930 , Heat transfer to liquid metals m 
concentric annuli, OM1-531 

Bender, I) J and Magee, P i T , 1969 , Heat transfer in 
rod Bundles with liquid metal coolants, GBAP-10052 

Biown, W J , 1968 , Heat transfer to liquid metals 
m double pipe heat exchangers, ANL-7 289 

Deissler, R G and Taylor, M F , 1957 , Analysis of 
axial turbulent flow and heat transfer through bank, 
of rods or tubes, PID-7529 (Pc 1), Book 2, 416 

Dwyer, 0 E and Tu, P S , I960 , Analytical study of 
heat transfer rates for parallel flow of liquid metals 
through tubes bundles, J Heat Transfer, 86 


Friedland, A J and Bonilla, C F , 1961 , Analytical 
Study of heat transfer rates for parallel flow of 
liquid metals through tube bundles, A I Ch ^ J f 7 


Lyon, R H (ed\, 1955 . Liquid Metals Hand Boolg, US 
Atomic Energy Commission 


Mill saps , K and Pohlhausen, K. , 1956 . Proceedings of 
the conference on differencial equations, Umv of 
Maryland, Maryland 

Nijsmg, R , and Eifler, W , 1973 , Temperature fields 
in Liquid metal cooled rod assemblies, Progress m 
Heat and Mass Transfer, Vol 7 


Ober John, ¥ J , 1970 , Turbulent flow thermal -hydraulic 
characteristics of hexagonal pitch fuel assemblies, 
Liquid metal heat transfer and fluid dynamics, ASME 


Raznjewic, H , 1976 . Hand hook of thermodynamic tables 
and charts, McG-raw Hill 


Schneider, P J , 19 56 , Heat transfer fluid mech 
Stanford Umv Press, Stanford, California 


f 


12 


Inst 


61 


13 Singh, b H , 1957 » Appl Sci Res A 7 

*^■4 oleicher, 0 A and Tnbus H , 1957 » Heat transfer 

m a pipe with turbulent flow and arbitrary wall 
temperature distribution , ASMS 79 

15 Stem, R P , 1966 , Advances m heat transfer, Vol 3 

16 Rays, W M , 1975 . Convective Heat and Mass Transfer. 
Tata McG-raw Hill 

II Ye vick, J G , 1966 , Past Reactor Technology Plant 

Design, MIT Press 


l 


/ 

4 


** 

4 




appel di a a 


[A 1] SEPARATION OP VARIABLES 

The solution of Eq (3 12) and (3 13) involves H, K, 
Ky and L as nondimeneional parameters The solution 
requires the knowledge of f 1 ( r 1 ) an< ^ S 1 ( r 1 ) 13 

determined both by Prandtl and duct Reynolds number Values 
of these nondimensional parameters determine for a specific 
heat exchanger, what may be called as operating conditions 
for gheat exchanger 

To solve Eq (3 12) and (3 13) separation of variable 
method is adopted Let, 

T.Jr^ z) = B^r.^ ©(z) 

substituting Eq (A 1) in the Eq (3»12) and (3«13) results 
m Eq. (A 2) and Eq (A3) 


Tube 


1 

^ . A. 

©( z) 


a©(z) 

a 

dz 





“i&> 

dr. 


(A 2) 


Annulus 



d©( z) ^ 

ii ” g 2 (r 2 ) “ 2 B 2 Cr 2 )(i 2 +a) ii 2 


* dE 2 (r 2 ) 

[(r 2 + a) f 2 (r 2 ) - 2 ~ ] 


dr. 


(A 3) 



N> 


A- 2 


The left hand sides of the above equations are functions of 

A 

only whereas the right hard sides are functions of r^ 

only , consequently, they can be equal only if they each 

2 

equal a constant, - X 


d©( z) 


- x' 


(A 4) 


e(z) d(z) 


2 * 

-X z 

or e(z) = e 

The right hand side of Eq (A 2) and (A 3) can 


(A 5) 


dr. 


[*!<*!> ? 1 


dE, (r, ) 2 * * * * v n 

i + * 6i( r i) r i B i ( r i' = 0 
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( A. 6) 


and 


A dE (r ? ) 

— [f 2 (r 2 )(r 2 +cr) x ] 

dr ? 2 2 2 dr 2 


- aw 2 g 2 (r 2 ) (r 2 +o)E 2 (r 2 ) 

(A 7) 


The interior boundary conditions become, 


< 11 ^( 0 ) 


= 0 


dr. 
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dE 1 (l) 
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dE 2 (0) 


dr. 




dB 1 (l) 

dil 


S 2 (0) - E 1 (l) 


(A.8) 

(A.9) 
(A 10) 



0 


A-3 
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dB (1) 



(A 11) 


[A 2] TWO BA&IOB bTUBIvi-IIOUVILlJi PROBLEM 

The Eq (/ 6) and (A 7) correspond to what may he 
called a two region Sturm-Liouville problem and may be 
written as a single differential equation 

sl (k dk'* + ^ & = 0 » 0 lx <_ 2 (A 12) 

with boundary conditions, 

y' (0) = 0 and y' (2) = 0 (/« 13) 

where x = 0 and x = 2 correspond to r-^ = 0 and r 2 = 1 
respectively 

The functions y, k, g are defined over both regions, 
tube and annulus, of the heat exchanger by, 

Tube 


(x 

A 

* r l 

f o < r x < 1 ) 


y 

* 

= E x ( 

:?i) 

(A 14) 

k 

A 

= r l 

■^i ( r i ) 

(A 15) 

S 

A 

= r l 

gi ( ^i ) 

( A 16 ) 
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Annulus C - 

r 2 t 1 , 

0 1 ? 2 5 


y = 

E 2 (r 2 ) 


(A 17) 

h 

(r 2 + a) 

f 2 ^ r 2 ^ 

(A 18) 

g 

- to (r 2 

+ o) g 2 (r 2 ) 

(A 19) 


where the function y must meet tne compatibility conditions 
at the inner wall of heat exchanger Eq ( A 9 ) and (A 10) 
and the functions k and g have finite discontinuities at 
the inner wall of heat exchanger 

The difficulty m applying the Sturm-liouville theory 
to the present system lies m the fact that it specifies 
that the functions must be continuous The functions are 
not continuous m the system being considered Because of 
the fact that the functions are continuous m their regions, 
1 e , piece-wise continuous, Stem (1966) has been able 
to show the Sturm-Liouville theory applies to the present 
system 

According to Stem, then, there are infinite sets of 
eigenvalues and corresponding eigenfunctions, having limit 
points from - «* to ® Thus m one region problem the 
negative eigenvalues do not exist but for two region 
problem, having discontinuity at interface, both negative 
and positive eigenvalues exist 



The E'q (A 4), (A 6) and (A,7) according to above two 
region Sturm-Liouville problem have infinite sets of positive 
and negative eigenvalues and corresponding eigenfunctions 
These equations thus may be written as, 


_ , 2 
e(z) dz “ “ 11 

A or e to » t 

d _ / A v A V, ( r T ) 


(A 20) 
t A Zi) 


* 1 J + X n 2 *L ^A) E l,n<*l> = 0 


(A 22) 


-4- [i 2 (A.) (V»> dB?l ^ (r2) 

dr 2 dr 2 


- * n 2 « 2 (J 2 +o) g 2 (J- 2 ) 5 2>n (i 2 ) =0 (A 23 


and boundary conditions 


dE (0) 


(A 24) 


dE l,n<°> A dI W°> 


(A 25) 


“l.n^ 


E 2,n<°> - E l,n^ 


(A 26) 


dS 2.n fjj 


(A. 27) 


n = + 0, + 1, + 2, , +ob 


I 


A -6 


«4» 

where is the eigenfunction corresponding to X^ 

[A3] general SOLUTION 

The general solution to the Eq (A 6 ) and (A 7 ) may 
he written for region 1 and 2 as follows 


Tube 


A A 


n = + « 


. 2 - 


11 (r l’ 2) = a J °n h.n <*!> • " 


(A 28) 


Annulus 




n I!I 5 » "*2.» (r V • 


2 A 
-*n 2 


(A 29) 


where Q n = Expansion coefficient associated with n-th 
eigenvalue 


In order to determine eigenvalues, eigenfunctions and 
expansion coefficients we first of all derive orthogonality 
property 

[A 4] ORTHOGONALITY OP EIGENPUNCIIONS 

The orthogonality property of eigenfunctions may be 
derived in the following way Bq (A. 22) is first multiplied 

A 

by E, „ to obtain a relation between different order 
J* p in 

eigenfunctions A second relation is obtained by inter- 
changing m and n m first relation. The two expressions 
are then integrated between 0 and 1, subtracted and 
simplified to yield Eq (A 30) below 


A-7 


(X 


n * X .*> B l l V <*!» 


= 2E- 


dE 1 (1) * dE, (1) 

(1) [ -i^S B l, m CD — ^ 1 U 30) 

dr, dr. 


A similar treatment to Eq (A 23) yields Bq (A 31) given 
below 



x * 2 > m 2 {S 2l „( ? 2 ) 

2 r A ^2 yJO) * (0) 

= TT C L 2,m (0 ^ ~-f* E 2,n (0 > ~ ^ 3 U 51) 


where n = + 0, + 1, + 2, 

The Eq (A* 30) and (A 31) are combined by applying the 
innerwall boundary conditions to yield 


A A 


( X n ~ X m ) [B x tSq^C^) 


- HB 2 d 2tn { r 2 ) B 2>m (r A 2 ))] 0 


(A 32) 


m 

= ± 0, + 1, + 2, 

n 

This gives desired orthogonality condition 


B 1 { \,n< r l> - ®2 (B 2 >n < ? 2 ) E 2 >t A)> 3 °> 

n ^ m ( A 33) 
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B l ([S l,n (# l^ 2> - EB 2 {[S-fh] 2 ) = 


2n v ' t '2 J ' = N n » n = m 


U 34) 


n 


+ 0) + 1, + 2 , 

where { } and ^2 { } physically mean average over cross- 
section They are defined as follows 

Tube 


B i { > = 2 / g 1 (r 1 ){ } r ± d^ 


(A 35) 


Annulus 


B 2 { } 


1 + 2a o 


Mis 

I { } (r 2 + cr) d *2 (A * 56 > 


Thus the average temperature over the cross-section (bulk 
average) is given as { T i (r^, z ) ) * which is denoted as 

T iav ^ z ) 

Note that dimensional equivalent for hulk averages 
are defined for tube and annulus, as follows 


Tube 


V > 


;i2 12 

= (2t ij u i( r i) r i { >dr 1 )/(27t |u,Ct)r, dr^ ) 

0 0 


(A 31) 
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Annulus 

r 22 r 22 

B 2 { } = ( 2k J u 2 (r 2 )(r 2 ) { > dr 2 )/( 2 k Ju t L^r dr ) 

r 21 r 21 2 

( A 38 ) 

[A 5] EIGENFUNCTIONS 


Solution of the Eq (A 22) and (/ 23) for \ Q 2 = 0 

are E l,0 anE E 2,0 e< 3hal to a constant No generality 

is lost m assuming those constant equal unity 


A A 



The zeroth order normalizing factor as given by Eq 
then becomes, 


(A 39) 
(3 61) 


N 


o 


= 1-H 


( ' 40) 


solution of Eq (A 22) and )(l 23) is best effected by 
defining auxiliary functions F(r 1 , x^) and G(r 2> X n ) Let 




li 


2,n^ t 2 ) 





where F(r 1( X^) and G(r 2 , X n ) are solutions to 


U 41) 
U 42) 


dr. 




dF(r, 


V 


dr. 


] + X 


n 


Sj.(r 


l) y(r L , 




(A 43) 
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and, 


tf 2 (f 2 ) (? 2 +°) — 1 ~ ■— J - 

ar o dr„ 


2 . 2 


A A 


X n w s 2^ r 2^ r 2 +a ^ S ( r 2 ? x n^ ~° ( A 44 ) 


In;Liial..CrQiidj. tio ns 


di(0, X n ) dF(l, Xj 

- = 0 j ; 

dr n dr„ 


= 0 


F(0, x n ) = 1) G(l, Xj = 1 , 


n' 


( 1 45) 
(A 46) 


The constants are so defined that they satisfy inner vail 
boundary condition Bq (A 25), 


K 


dE., (1) 

l,n v ; 

A 

dr n 


dE 2.n (°) 

A 

dr. 


(A 25) 


1 “*2 
which in terms of auxiliary functions become, 


KA* 


dP(l, X) 


n 


n 


B, 


, aG(o f x n ) 


dr. 


n 


47) 


dr. 


Therefore let, 



< 

aa(o, x n ) 
a? 2 

(A 48) 

and 

! 

B 

n 

t* 1 ’ v> 

* 

(A 49) 


dr 2 


where 

n - 

Hh 1, + 2, + 3f 

(A 50) 
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[A 6] EIGENVALUE EQUATION 


iho eigenvalues are determined by the innerwall 
boundary condition Eq (A 26) which m terms of auxiliary 
functions P and G becomes 




ag(o, Q am, x n ) 


di . 


A 

dr. 


Iv 


«(i, K) 


n' 


dr n 


G(0, K) 


n' 


dG(0, X ) 

— P(l, X n ) (A 51) 


A 

dr. 


n = ± !, i 2, + 3, 


The eigenvalue equation is therefore defined as, 



dG(0, X ) dP(l, X ) dG(0, x) 

\ — t» ^r- 2 - + — » n ) 

drg dr^ dr g 


The eigenvalues 
equation, 



dF(l, xj 



G(0, X n ) 


U 52) 


are obtained fxom the roots of the 


U ) - 0 

' y n' 


(/. 55 ) 


with 


0 


A-12 


[ A 7] hXSlM SION COEFFICIENTS 

Temperature distributions m tube and annulus are 
given by Eq (A 28) and (A 29) 


Tube 


*i z) 


0 1 (A) e "" ■ (' a) 

l n l.n 1 
n= -« 


. 2* 
-X. z 


A nnulug 


A A A 

T 2 (r 2 , z) 


n= +« /. * * 

y C (r 2 ) e 

h i.» n 2,n 2 


, 2 A 
' X n 2 


(i 29) 


Relation defining expansion coefficients, C n , are obtained 
from the entrance and exit temperature profiles, of the heat 

exchanger fluids 


z = 0 

T l( r l' 
T 2 (r 2 , 0) 

Z = L 

*1<V 

*• a * 

T 2 (r 2 , 1) 


0 

P 2 (r 2 ) (unknown) 

P 1 (r 1 ) (unknown) 
1 


(A 54) 
U 55) 


(A 56) 
(A 57) 
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The end temperature distributions may be written m terms 
of the series solution, 

A 

z = 0 


+ ‘ 


0 


n=-" n E]L » n ^ 1 ‘ 


(A 58) 


-f. 00 


= n | Jn t ‘2,n < r 2> 


(A 59) 


z 


±* - 


2 

- X I 


Pl( r l> - I Onh.n (r l> e “ 

n» - « 


U 60) 


-f. 00 


l C E_ (r 2 ) e 
n=-» 11 2 ’ n 2 


2 * 
- X ^ L 
n 


(A 61) 


Manipulation of above equations yields expressions for the 
expansion coefficients Eq {n 58) and {a. 59) are multiplied 
by EL and HB 0 . respectively, and the bulk average operation 
is applied to the result The expressions are then subtracted 
and the orthogonality condition, Bq (A 35) and (A 34) applied 
to give Bq, (A 62) below A similar treatment to Eg ( 60) 

and (A 61) yield Eq (A 63) 


- - m 2 tB 2,n (? 2> f2 (r 2 >} 


( A. 6 2 ) 


" n 

0 N e n 
n n 


. 2 C 

X„ L 


A A 




-HB 2 (B 2in (r g )> 


(A. 63) 


H 



A-14 


the above expressions cannot "be used directly as they 

(r 2 ) However, 
these unknown functions can be eliminated by substituting 
them as their series expansions After substitution the 
two Bq (A 62) and (A 63) are added to give Eq (A 64) 


contain the unknown expressions p-^r^) and P 2 


A 



N. 


n 


(1 + 




x B 1 { S.,n ( r l) ®i,k ^ r l^ } " m 2 ( ^2,n^ r 2^2,k ( ' r 2^ } 


- HB 2 (B 2 >n (r 2 )> 


(A 64) 


n 


0, + 1, + 2, ± 3, 


For n 4 0, tho term for K = 0 drops out of above summation 
and 0 Q does not appear m the expression An explicit relation 
for m terms of higher-ordered expansion coefficients may 
be obtained for n = 0 


Foi n e 0, Eq (A.64) becomes, 


A 



N. 


k = -» 


k 


(e 


2 *- 


B 1 {E l,k (r l )} 


- H 

(A 65) 


A fairly straightforward manipulation of the Bq (A 64) 
yields the following simplified form, 


i 
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4 

M 

f 


I Gjc ( l~e 
k= 


-x 


2 

1a 


I 

) 



lr 



A 



< ; i>> 


(A 66 ) 


where, 


*a,k “ B 1 <B l,n ®i, k <*i)> > * 4 n (i 67) 


« 


n.n = [V£V (r l>J 2,+ < e 


2 A 


1-e 


p A 

2 i 


■) N n ] , k = n 


(i 68) 


n = + 1, + 2, + 3, 


The zeroth order expansion coefficient does not appear in 
the above equation and must be calculated by Eq (A*65 ) 

[A 8] NORMALIZE! EIGENEUHO TI0N8 


At this point it becomes desirable from practical 
viewpoint to normalize the eigenfunction with respect to the 
square root of the absolute value of N n Absolute value is 
used because for n < 0, N n becomes negative This normaliza- 
tion is desirable for computational considerations, since 
it was found that unnormalized eigenfunctions become 
unmanageably large as their order is increased let 



(A 69) 



S 2,n ^ 2 ) 
*■»' ■* . ■■■■—■■ ■■ ■ ■ 


/ l I nl 


(A 70) 
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\11 the eigenfunctions are normalized except zeroth— order 
eigenfunctions, which remain, 

E l,0 = E 2, 0 = 1 

using normalized eigenfunctions, the orthogonality condi- 
tions become 


B i tL i, n <n> - ™ 2 E 2 im ( ? 2 )> = 0 > 


B. 


n 4 m 
2 


1 <[%,n - HB 2 {CB 2>n (r 2 )]^> = 


U 71) 
+1 n > 0 
-1 n < 0 


n 


( A 72) 

Relations which were derived so for one still valid with the 
value of H n taken as + 1 for n > 0 and - 1 for n < 0 

[A 9] DEFINITIONS TO AVOID USE OF NEGATIVE INDICES 


Eigenvalues are given as 




x 


2 

m-1 ’ 


X 


2 

m 


—o> < m ^ “ 

MW 1— ■ 

2 

Wo define auxiliary eigenvalues P n 


where 


Such that 


X 2 -v - p 2 
A ~n p n 


(/ 73 ) 


Eigenfunctions associated with negative indices are 
also redefined as, 




— n (^) = 


1,-n 


E l,n (-l) 


U 74) 


B 2,i-n ^ r 2^ 


E 2,n ^ r 2^ 


(A 75) 


Expansion coefficients associated with the negative eigenvalues 
may be redefined in terms of auxiliary expansion coefficients 


and a n 


0 + c 

n n 

2 * 

-e 1 

- <• n >*n 


(A 76) 
( * 77) 


2 2 

The above definitions now give summation for X n and p n 
from 0 to +“ and thus obviates use of negative indices 


APPENDIX B 


UNIFORM HEAT FlUX SAMPLE CALCULATION FOR 
XALPAXXAM BATA 


[B 1] NUMBER OF HEAP TRANSFER UNITS 

The number of heat transfer units (NTU) for a 
heat exchanger is defined as UA/C min where UA for double 
pipe is defined from equation (B 1), 


1 

UA 


2flk L 
w 


2 \ 


2k. 


'w 


h^Ag 


+ in (^)J 

r 22 


(B 1) 


and G = Minimum heat capacity flow rate 
mm 

A^ s Surface area of tube 
A 2 = Inner surface area of annulus 
zs Conductivity of tube metal 

- Heat transfer coefficient for tube side 

- Heat transfer coefficient for annulus side 
= Length of heat exchanger 


h l 

h 2 

L 


For double pipe heat exchanger the following 
convenient expression for HUT, in terms of Nusselt numbers 
can be derived from Eq (B 1), 




c , 

min 

2*V 


2k. 


w 


Nu 1 k 1 


3w ( r ?2~ r P1 ) 
Nu 2 k 2 r 21 



(B 2) 


The Fusselt numbers were obtained using experimental 
correlations of Loyn (1955) for tube and Bailey (1950) for 
annulus 


Tube 

Hu l = 7 + 0 025 P ei ° 8 (B 3 ) 

Annulus 

Nu 2 a 5 25 + 0 0188 Pe 0 8 ( — )° 5 (B 4) 

r 22 

In Bq (B 2) various parameters for JCalpakkam heat exchanger 
have following values 


Ku 1 a 8 44 

Nu 2 c. 6 676 

\ = 19 77 W/m°C (at 425°G) 

= k 2 =» 80 watts/m°G )(at 425°C) 

G min * 110 watts 
L * 2 48 m 

r 2 ^_ = 11 x 10”5 m 

r 12 ~ 7 x m 

and 2 ( r 22** r 21^ a H y draulic diameter 

4 X Gross -sectional area of Annulus 

) 


( 


Wetted Perimeter 


B 3 


With these yalues Eq (B 2) gives 

WTU = 10 470 

[ B 2] EFFICIENCY 


Efficiency for a heat exchanger is defined as the 
ratio of actual heat transferred uo that which would "be 
transferred if heat exchanger were infinitely long for 

counter flow heat exchanger, efficiency in terms of NTU 

nnd ft fu _ P rl -ary Heat Capacity Fl ow Rate 

Secondary heat Capacity Plow Sate ^ 1S g;rv ’ en as 


H < 1 

1 - Exp (-NTU (1 ~ H)) 

e = (B 5) 

1 - H Exp (-NTU ( 1-H) ) 

H > 1 

1 - Exp (- NTU (1 - -rr)) 

e £ (B 6) 

1 - Exp (-NTU (1 - —)) 

Prom equation (B 6) efficiency can be calculated substi- 
tuting the values of NTU and H For present case H = 1 65 
Therefore Eq (B 6) gives 


e 


0 993 


[B 3J OUTLET TEMPiifiATUkt, 8 


B h 


1 


i 

1 



Efficiency for a given heat exchanger m terms of 
temperatures is defined as follows 

H > 1 


e ~ 


t il ~ T 10 

T 2L ~ T 10 


(B 7) 


H < 1 


T 21 T 20 

e = T ZTv — 

21 1 10 


(B 8) 


where subscript 

1 denotes tube side and 2 denotes annulus side, 

0 denotes entrance side and L denotes exit side 

From general heat balance we can write, 

(T 1X, - T 10> = H < T 21 - T 20> < B 9) 


For present case secondary outlet temperature, T^, can 
be determined from Sq (B 7) liimary outlet temperature, 
T 20> can be determined from Eq (B 9), They are given 
as , 



518 0 0 


377°0 


/ 


i/ 
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